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Abstract 

We demonstrate that in Optimal Velocity Model (OVM) delay times 
of vehicles coming from the dynamical equation of motion of OVM almost 
explain the order of delay times observed in actual traffic flows without in- 
troducing explicit delay times. Delay times in various cases are estimated: 
the case of a leader vehicle and its follower, a queue of vehicles controlled 
by traffic lights and many-vehicle case of highway traffic flow. The re- 
markable result is that in most of the situation for which we can make a 
reasonable definition of a delay time, the obtained delay time is of order 1 
second. 
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1 Introduction 



The history of traffic dynamics began early in 1950s. Two different models were 
proposed: the one is car- following modelUT], 0, ^, which describe the motion 
of vehicles by many-variable differential equations, and the other is the fluid 
dynamical model [Q, which regard traffic flow something like fluid. These two 
models have been further investigated by many authors. However either model 
did not succeed in explaining the behavior of real traffic flow in some points. 
Especially those traditional models can not attain an unified understanding the 
most remarkable fact that traffic flow has two phases; one is free flow with low 
car-density and high average speed and the other is congested flow with high 
car-density and low average speed. 

Recently, in these two models, there are found new mechanisms which explain 
the existence of two phases. In car-following models, a modified model was pro- 
posed by introducing optimal velocity^. (We call this 'Optimal Velocity Model' 
(OVM) hereafter.) The equation of motion of OVM is given by 

x„(t) = a{y(Ax„(t))-i„(t)} n = l,2---Ar, (1) 

where the notations are; car number n, time t, position of n-th car Xn and its 
headway Ax„. A dot on the variable denotes differentiation with regard to time 
t, sensitivity a is a constant parameter. The essential difference from traditional 
car-following models is the introduction of an optimal velocity V{Ax) of a vehicle, 
which value is a function of headway distance. A driver reacts according to the 
difference between the vehicle's velocity and the optimal velocity V{Ax) and 
controls its velocity by accelerating (or decelerating) his vehicle proportional to 
this velocity difference. The dynamical equation of OVM has two different kind 
of solutions. One is a homogeneous flow solution and the other is a congested flow 
solution which consists of two distinct regions; congested regions (high density) 
and free regions (low density). In OVM, if the density of vehicles is above some 
critical value the traffic congestion occurs spontaneously, which can be understood 
as a sort of phase transition from homogeneous flow state to congested flow 
statel, |]. 

Also in fluid dynamical models, the following equation has been proposed]^. 

dv dv V{p) — V cldp P d'^v . , 

dt dx T p dx p dx"^ ' 

where Cq, / and r are constant parameters of the system and density p and velocity 
V of vehicles are functions of location x and time t. They call the function V{p) 
'safe velocity'. To see the similarity of these two models, let us ignore second and 
third term of rhs of Eq.(^ and use D/Dt = d/dt + vd/dx; the differentiation 
acting on the variables of individual vehicles. The fluid dynamical equation (|^) 
becomes 

Dv _ Vip)~v 

Dt ~ T • ^"^^ 
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This is quite similar to Eq.(|T|) with optimal velocity V{Ax) and sensitivity a 
replaced by 'safe velocity' V{p) and parameter 1/r, respectively. The reason why 
these two models can explain the formation of traffic congestion comes from the 
form of dynamical equation and the introduction of 'optimal' or 'safe' velocity. 

Contrary to this, the original equation of motion of traditional car-following 
models PI 

Xn{t + r) = X{Xn-l{t) - Xr,{t)} , (4) 

is essentially a first order differential equation without a delay time r of response. 
Because traffic flow governed by a flrst order differential equation is always stable, 
the delay time r plays a crucial role in describing the behavior of traffic flow. 
The origin of the delay time r have been thought to be a physiological delay 
of response. In fact, it is well-known that the motion of a vehicle accompanies 
some delay time in response to the motion of its preceding vehicle. There seem, 
however, questions to be discussed: what is the role of the delay time r, or, 
whether or not the delay time r is equal to the value of observed one. 

Here we should remark that it is necessary to distinguish two different types 
of deflnition concerning the notion of "delay time" . There is a time lag until the 
driver begins an action after being conscious of a stimulus. It may also takes a 
flnite time for a vehicle to change its velocity after the operation of driver. We 
deflne such physiological and mechanical time lag as "delay time of response". 
This can be measured in principle although it may vary depending on the type of 
the stimulus, individual driver and performance of vehicle. On the other hand, 
we can deflne "delay time of vehicle motion" as a period from the time of velocity 
change of a vehicle to that of following one. The delay time of vehicle motion can 
also be estimated in observations of real traffic and/or numerical simulations using 
traffic models. However in general case, it is difficult to find a good deflnition of 
the delay time of motion of two successive vehicles. Only in the restricted case 
in which two vehicles behave quite similarly, for example, 

Vnit)^Vn+lit-T), (5) 

we can deflne the delay time of motion as T. 

In the car-following models, there are another two types of "delay time" to 
be distinguished. One is the delay time r explicitly introduced as a parameter 
in the equation of motion (see Eq.(^)), which we call "explicit delay time" in 
this paper. This may correspond to the delay time of response. The other is the 
delay time emerging as a result of the dynamics of traffic flow, which is quite 
different notion from the explicit delay time. This will correspond to the delay 
time of vehicle motion stated above. However, this delay time includes not only 
the contribution from the explicit delay time but also that from purely dynamical 
origin. We call the latter as "effective delay time" . Obviously, the effective delay 
time is not zero even if we introduce no explicit delay time. 
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In this paper, we discuss the delay time of motion in OVM with no exphcit 
delay time. Therefore, because the resulting delay time is equal to the effective 
delay time, we can investigate the delay time from dynamical origin. This quan- 
tity may be compared with observed delay times of vehicle motion in various 
cases. It will be clear that the effective delay time obtained by our procedure is 
enough to explain the delay times observed in actual situations of traffic flow. 

In section 2, we define the effective delay time of vehicle motion in OVM in 
terms of the vehicle motions of a leader and its follower and make an analytical 
study within linear approximation. Then we carry out numerical simulations to 
obtain the effective delay time in several cases. We show the results in actual 
traffic situations: the effective delay times under traffic rights (section 3) and 
those in uniform traffic flow and in congested flow (section 4). Discussions and 
further prospects are given in section 5. 

2 Delay Time for Leader and Follower Case 

First let us make an definition of effective delay time of vehicle motion. Consider 
a pair of vehicles, a leader and its follower. This pair of vehicles may be either 
separate two vehicles or any pair of vehicles picked up from a series of vehicles in 
highways or from a queue waiting to start under traffic lights. 

When a leader moves with its velocity v{t), and its follower replicates the 
motion of the leader with some delay time T (then the follower's velocity is given 
by f (t — T)), we can define the delay time of vehicle- motion as T. It must be 
remarked that we do not define the delay time as the delay of motion of the 
follower with respect to its position but its velocity replication. 

Let the positions of a leader and its follower be y{t) and x{t). In this case 
Eq.(|l|) is written as 

xit)=a{V{yit)-xit))-x{t)}. (6) 
Uniform motion is described by 

yo{t) = V{b)t + b, xo{t) = V{b)t, (7) 

where b is headway and V{b) is a constant velocity. To investigate the response 
of the follower vehicle to the leader vehicle, we introduce a small perturbation 
\(t) and its response ^(t): 

y{t) = yoit) + A(t), x{t) = xoit) + iit). (8) 

Inserting Eq.(|D into Eq.(|^) and taking a linear approximation, we get 

m+am+afm = cifm. (9) 

where / = V'{b). This is just equivalent to well-known equation of motion for 
forced oscillation with damping term caused by friction. 
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In order to find a solution, we first write X{t) by a Fourier expansion 

X{t) = J A(cj)e*^*dcj . (10) 
For the Fourier component Aoe*"^*, the solution of Eq.(P) is given by 

m = , e^^*. (11) 



This is rewritten as 
where 



m = |r/|Aoe-(*-^), (12) 



I.P = '^'^ (^^)^ (13) 

XI (a/-o;T + M'' ^ ' 

T = itan-i-^. (14) 

When / < a/2, the amphtude |?7| is monotonically damping function of uj. On 
the other hand, when / > a/2, |?7| takes its maximum at u = Uo', 

ul = a{f-a/2), (15) 

so we call this ojq as "enhanced mode" (see Fig.0(a)). Note that / = a/2 is the 
critical point for the instability condition of homogeneous flow as we found in the 
previous paper [Q. Eq.(|T5D shows that the enhanced mode ojq 0) exists so far 
as the instability condition / > a/2 is satisfied. 

Let us examine some characteristic cases. For / < a/2 where low frequency 
modes dominate \uj\ <^ a, f, we have 

|?7| ~ 1, ^ ~ y- (16) 
In this case the response ^{t) to the perturbation A(t) becomes 

e(t) = J \{uj)e'^^'-^^duj = X{t - T), (17) 

which leads 

x{t) = V{b) + i{t) = V{b) + X{t -T) = y{t - T). (18) 

Thus for sufficiently slow perturbation, the delay time of vehicle motion becomes 
T in Eq.(|T6|), which is approximately the inverse of derivative of the optimal 
velocity function at corresponding headway. 
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In the other case / > a/2, the amphtude \t]\ takes maximum value at u = ujq. 
Then we have 

T — J_f -1 na^i 

-'enhanced — tan , {^^) 

ujq a 

which indicates that the delay time T for this enhanced mode depends on the 
sensitivity a in contrast to the previous case Eg . ([TBD . One can easily confirm that 
T tends to 1// when a is close to its critical value 2/. We should remark that 
an exact replication indicated in Eq.(|18D is not realized in this case because |?7| 
is not always equal to 1 in Eq. (]T^) . (See also Figs.|I] and the discussions below 
Eq.(ll).) 

Now let us see the results of numerical simulations and compare them with 
those of our analytical consideration. Here and hereafter we use the following 
form of the optimal velocity function P], 

V{Ax) = 16.8 [tanh 0.0860 (Ax - 25) + 0.913] (for Ax > 7m) (20) 
= (for Ax < 7m) (21) 

whose parameters are determined from the car-following experiment on Chuo 



Motorway [10, 11]: the inflection point is (Ax, x) = (25 m, 55 km/h), the maximal 
velocity is Knax = 115 km/h and the minimal headway is Axmin = 7.0 m, which 
includes the length of the vehicle (5 m) which was used in the experiment. 

In numerical simulations, we prepare a pair of vehicles with their unperturbed 
motions of Eq.(^ with headway b. If the leader changes its motion by X{t) then 
the function ^{t) can be obtained by the numerical simulation from which the 
delay time can be read off. We choose X{t) of the leaders motion as 

y{t) = V{b)t + b- Xo cos iut, y{t) = V{b) + Xqu sin iut for t > (22) 

with various u setting uXq = 0.1 m/s, which means that \y(t) — V{b)\ < 0.1 m/s. 

As illustrations, we show the behaviors of X{t) and its response ^(t) for a = 
2.0 s-\ b = 25 m (therefore Uo = 0.938 s"^). Fi gures are the cases for u = 
0.1, 0.938, 1.5 s~^. To find the values T, we first rescale ^{t) and then translate 
it so as to coincide the curve X{t) in Figures 0. The value of jr/l is this scale factor. 
The numerical simulations have performed for 6 = 25, 30, 35, 45 m. In Figures 
1^, we show the numerical values for \7]\ and T. The corresponding analytical 
results Eqs.(|l3D and (|14D are also shown in these figures. In both figures, the 
analytical and numerical results agree quite well. 



3 Delay Time under Traffic Lights 

The delay time of vehicle motion is clearly recognized in motions of a series of 
vehicles under traffic lights. Consider the situation in which every vehicle waits 



5 



15.45 



15.4 



15.35 



15.3 



15.25 



15.2 



15.15 




15.45 



15.4 



15.35 



15.3 



15.25 



15.2 



15.15 



10 20 30 40 50 60 70 80 90 100 
time (s) 




10 15 20 25 30 35 40 
time (s) 





15.4 p 




15.38 ■ 




15.36 - 




15.34 1 




15.32 1 


i 


15.3 






'G 
_o 


15.28 - 






> 


15.26 ■ 




15.24 ■ 




15.22 - 




15.2 - 




15.18 - 




Figure 1: Numerical results for the motion of leader and follower, where 
Ax =25 m and a = 2.0 s~^. Frequencies of leader's motion are (a) uj = 0.1 s~^, 
(b) UJ = uq and (c) u = 1.5 s~^. 



until a red light changes to green. The initial conditions are as follows; 

Axi(0) = oo, i;i(0) = 

Ax„(0) = 7(m), x„(0)=0 (n = 2,3,...) (23) 

If the light changes to green, the top vehicle will first accelerates to start, followed 
by the succeeding vehicles according to the equations of motion. 

We perform the numerical simulation for the cases a = 2.0, 2.8 s~^ and obtain 
the time dependence of velocity x„(t) of n-th vehicles in a queue. Figures H(a) 
and ^(a) show the behavior of motion for the first ten vehicles. In the figures, 
it seems that the behaviors of motion for vehicles on the down stream converge 
into a common shape. Thus, except first several vehicles, every vehicle in the 
queue almost replicates the behavior of its preceding one with a certain delay 
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Figure 2: Each curve shows the behavior of (a) \ri\ (Eq.([T3|)) and (b) T (Eq.(|T3|)) 
for h = 25, 30, 35, 45 m with a = 2.0 s~^. Plotted marks on the curves show 
numerical results. 

time. In this case, the relation can be applied for n-th vehicle with enough 
large number n. In this sense we may say that replicative pattern of vehicle 
motion is realized "asymptotically" and only in this occasion we can define a 
delay time of motion T in just the same way as we defined in section 2. Note 
that this definition of the delay time slightly differs from a time lag with which 
the successive vehicles start for example. 

From Figures |^(a) and §(a), we can estimate the delay time T = 1.10 s for 
a = 2.0 s~^ and T = 1.03 s for a = 2.8 s^^. To confirm this similarity of vehicle 
motion, we also show plots of translated data x„(t — {n— 10)T) for the 7th, 8th, 
9th and 10th vehicles in Figures ^(b) and ^(b). 

4 Delay Time in Highway Traffic Flow 

Next we investigate delay time under a simple situation where vehicles move 
on a single lane circuit with circumference L. Of course we assume that road 
conditions are uniform along the circuit and drivers are identical. Numerical 
calculations are made with the initial condition: x„ = 26ni—b{n—l) {1 < n < N), 
Xn = V{b)- The first term of the right hand side of this condition means that a 
small perturbation is added to the first vehicle xi. 

In the previous papers we have shown a homogeneous flow changes into 
congested flow spontaneously if the density of vehicles is greater than the critical 
value. The results of simulations indicate that after enough time the traffic flow 
on a circuit creates an alternating pattern of high and low density regions. The 
motion of vehicles in this flow is visualized by plotting them in the 'phase space' 
(Ax, x). After the traffic flow becomes stationary, the trajectory of every vehicle 
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time (s) time (s) 

Figure 3: (a) Behavior of velocities of first ten veliicles under traffic lights with 
a = 2.0 s~^. (b) Figure of shifted curves Xn{t — {n — 10)T) with T = 1.10 s. 




time (s) time (s) 

Figure 4: (a) Behavior of velocities of first ten vehicles under traffic lights with 
a = 2.8 s~^. (b) Figure of shifted curves Xnit — {n — 10)T) with T = 1.03 s. 

in this 'phase space' draws a kind of limit cycle which we named 'hysteresis loop' 
in Ref.0 (see also Fig. 0). 

Now let us estimate delay times for two cases under this traffic flow: (A) the 
first stage and (B) the final stationary-state stage. 
■ Case A 

In this case the traffic flow is almost homogeneous. Let us pick up a pair 
of vehicles n = 10, 11. A small perturbation of the first vehicle Xi propagates 
backward and after several seconds those pair of vehicles change their velocities. 
The typical behaviors are demonstrated in Figures ^. In the same way as section 
2, the delay time of motion can be estimated from numerical results. The obtained 
values of delay times are shown in Table |I]for the case of a = 2.0 s~^ and Table ^ 
for the case of a = 2.8 s~^. As references, the delay times for low frequency limit 
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50 seconds and (b) Ax =40 m for 100 seconds with a = 2.0 s ^. 



and for the enhanced mode uq are also shown in Tables |l], ^ From these results, 
the low frequency limit is a good approximation and the delay time T is almost 
independent on sensitivity a for stable traffic flow. For unstable case, there exist 
contributions from blow-up modes and they have sensitivity dependence. 
■ Case B 

After sufficiently long time, traffic flow forms stationary patterns of high and 
low density regions. Under this situation vehicle does not change its velocity 
unless it encounters a boundary of high and low density regions. Let us observe 
the motion of a vehicle on the boundary. A vehicle which encounters the boundary 
changes its velocity. After some delay time, the following vehicle comes to the 
boundary and changes its velocity in the same way as the previous vehicle. The 
typical behavior of vehicles is shown in Figure ^ 

The delay time of vehicle motion in this case can be derived as follows. Con- 
sider two vehicles: one enters into a congested region from a free region and after 
a certain interval T, the next one follows. In the 'phase space' [Ax, x) (Fig.|^), 
the free region is denoted as a point F{Axp,vp), that is, vehicles is moving with 
velocity vp and headway Axp- The delay time T is deflned as the time which 
is needed for the next vehicle to reach the boundary and enter into a congested 
region. Now at the time when the first vehicle reaches the boundary of the free 
region the distance between the next vehicle and this boundary is of course Axp- 
The next vehicle runs with velocity vp and the boundary itself also is moving 
backward with velocity vb, so if the vehicle and the boundary will meet after 
time interval T, we can write the following relation 

vpT + vbT = Axp. (24) 

Similar relation can be written at a boundary where a vehicle exits from congested 
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Table 1: Delay times for various headway with a = 2.0 s^^. The second column 
indicates that traffic flow is stable (-) or unstable (+). The third and fourth 
columns show analytical results given in section 2. 
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Table 2: Delay times for various headway with a = 2.8 s^^. The second column 
indicates that traffic flow is stable (-) or unstable (+). The third and fourth 
columns show analytical results given in section 2. 
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region. 

vcT + vbT = Axc . (25) 

This can be confirmed if one recalls that the pattern of the flow is already station- 
ary, and the input vehicles of a boundary of congested region must be equal to 
the output of another boundary. The above two equations, (p4D and (^) indicate 
that the time interval T and vb can be graphically expressed as the slope and the 
intercept with vertical axis of the line connecting F and C for a = 2.0 s~^ and 
F' and C for a = 2.8 s~^ in Figure |^. If we write the velocity of the flrst vehicle 
as v{t), then the one of the following vehicle is v{t — T), which implies that T is 
just the delay time of vehicle-motion deflned in the previous section. 

It may be convenient to adopt the coordinate moving with the congestion 
pattern. Let x{t) and X{t) be the positions of a vehicle measured in a flxed and 
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Figure 6: Motion of successive two vehicles in congested flow. Initial condition is 
Ax =25 m with a = 2.0 s^^. 

a moving (with constant velocity vb) coordinates, respectively; 

x{t) = X{t) +VBt. (26) 

Then we get 



X(t) = a{ViAX{t))-{Xit)+VB)} 
= a{W{AX{t))-X{t)}, 



(27) 



where W{Ax) is the optimal velocity function in the co-moving frame; 

W{Ax) = V{Ax) - Vb. (28) 

Let us concentrate on the motion of vehicles in this co-moving coordinate. 
Every vehicle passes a same position with a certain time intervals T; in high 
density region, the time interval is given by 



Tr 



/AX 



c 



and (similarly) on the point F, we have, 

V X 



(29) 



(30) 



These two time intervals should be equal since otherwise the congestion pattern 
moves. We write this time interval by T; 



T = Tc = Tf 

which is of course equivalent to Eqs.(P^ and ( 
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Figure 7: The limit cycles for a = 2.0 and a = 2.8 s~^. 

We have carried out numerical simulations. For sensitivity a = 2.0 s~^, we 
find C(12.51, 2.05) and F(37.50, 28.55) yielding the delay of vehicle motion T = 
0.943 s and the back velocity of the boundary vb = 11-2 m/s. As for a = 2.8 s~^, 
C"(21.89, 10.92) and F'(28.11, 19.68) yielding T = 0.711 s and vb = 19.9 m/s. 

It is interesting to find that the main contribution of the resultant delay of 
vehicle-motion comes from the structure of the Optimal Velocity Model and not 
from the explicit delay r. 

5 Summary and Discussions 

The notion of delay time of response r has played a significant role in the history 
of traffic dynamics. Indeed delays of vehicle motions are observed in many cases, 
in traffic lights waiting queues or in highway traffic motions and the delay time 
are usually observed to be of order 1 second. We should take account of the effect 
of observed delay time and it has long been thought that it must be introduced 
in the equation of motion as an explicit delay time, most of which is caused by 
driver's physiological delay time and mechanical delay of response of vehicles. 
However it is known that the physiological response time is of order 0.1 second, 
not of order 1 second. We should be careful that the delay time of vehicle motion 
comes from another origin, that is, from the equation of motion itself which we 
have here investigated intensively. The results are summarized as follows: 

1. The case of a leader vehicle and its follower 

As is seen in Figure 0(b), if the headway distance is around 25 m in which 
drivers are sensitive to the behavior of the motion of the preceding vehicle, 
we clearly recognize delay time is around 1 second independently of the 
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frequency of the leader's velocity-change function A(t). However if their 
headway distance is more than 40 m, delay time is estimated to be larger 
than 1 second, and sometimes we obtain 6 second for the case Ax = 6 m 
for the case of low frequency (a; ~ 0). This is because of the structure of 
optical velocity function. If the slope of the optimal velocity function is 
very small, drivers are insensitive to the behavior of the preceding vehicle. 
This can be easily understood if one considers the extreme case in which 
the function V is independent of Ax (and so / = 0). In this follower 
never reacts to its previous vehicle and accordingly its delay time of motion 
becomes infinity. 

2. A queue of vehicles controlled by traffic lights 

In this case except the first several vehicles, most of the succeeding vehicles 
behave almost similarly as seen in Figures ^ and |. From those figures, 
delay times are read off; T = 1.10 s for a = 2.0 s~^ and T = 1.03 s for 
a = 2.0 s~^ Although T depends on the sensitivity adopted, the results 
obtained is again of order 1 second for reasonable realistic sensitivity. 

3. Many- vehicle case of highway traffic fiow 

In Figures § we show the typical behaviors of a pair of vehicles and in Table 
and 1^ delay times obtained by numerical simulations are summarized with 
various values of its headway Ax. Again in the case of Ax = 25 m, the 
estimated delay time is found to be of order 1 second, and for the larger 
Ax, the larger the value of delay time is obtained. In the case where the 
congested fiow becomes stable, T = 0.71 s for a = 2.8 s~^ and T = 0.9471 s 
for a = 2.0 s~^. Since the end point C (F) on the limit cycle becomes larger 
(smaller) as a becomes larger (see for example Fig. 7 in ref.|P]), so of course 
T becomes smaller for larger a (high sensitivity). 

All our results show that the estimated delay time in our OVM is almost 
enough to reproduce the order of observed delay time. It now becomes obvious 
that the delay time of motion arises as an effect of dynamical equation itself with- 
out any explicit introduction of r. It is interesting to find the main contribution 
of the resultant delay of vehicle motion are of same order of magnitudes in any 
cases. This may come from the structure of optical velocity function itself which 
we have determined phenomenologically. However we believe that this remark- 
able fact has its more profound reason, which will be made clearer by further 
investigation on structure of OVM by performing analytical study. 
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